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ABSTRACT 

Large-scale analyses of stellar samples comprised of cool, solar-like oscillators now com¬ 
monly utilize the so-called asteroseismic scaling relations to estimate fundamental stellar 
properties. In this paper we present a test of the scaling relation for the global asteroseismic 
parameter Vmax, the frequency at which a solar-like oscillator presents its strongest observed 
pulsation amplitude. The classic relation assumes that this characteristic frequency scales with 
a particular combination of surface gravity and effective temperature that also describes the 
dependence of the cut-off frequency for acoustic waves in an isothermal atmosphere, i.e., 
Vmax oscillations of cool main-sequence and sub-giant stars 

adhere to this relation, using a sample of asteroseismic targets observed by the NASA Kepler 
Mission. Our results, which come from a grid-based analysis, rule out departures from the 
classic g Aff scaling dependence at the level of ^ 1.5 percent over the full ^ 1560 K range 
in Teff that we tested. There is some uncertainty over the absolute calibration of the scaling. 
However, any variation with Tefr is evidently small, with limits similar to those above. 
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1 INTRODUCTION 

Cool main-sequence and sub-giant (i.e., solar-type) stars display 
solar-like oscillations, pulsations that are stochastically excited and 
intrinsically damped by near-surface convection. The solar-like ex¬ 
citation mechanism gives rise to a rich spectrum of detectable os¬ 
cillations. Data on global (or average) asteroseismic parameters as¬ 
sociated with the observed, rich oscillations spectra provide impor¬ 
tant diagnostics of fundamental stellar properties. Extracting these 
global parameters from oscillation spectra is usually very straight¬ 
forward, meaning a large number of stars can easily be analysed. 
This can be particularly advantagous when low S/N ratios either 
obstruct or prevent a more detailed analysis of individual modes in 
the spectrum. Moreover, detailed analysis of stars with high S/N or 
complicated oscillation spectra can be very time consuming. 

Two global asteroseismic parameters have received particular 
attention, and wide use. One parameter is the average large fre¬ 
quency separation, Av. It is an average of the observed frequency 
spacings between consecutive overtones n of the same angular 
(spherical) degree, /. The average large separation scales to a very 
good approximation as where p oc MjK is the mean density of 


a star of mass M and surface radius R (e.g. see Tassoul 1980, Ulrich 
1986, Christensen-Dalsgaard 1993). This scaling has physical justi¬ 
fication (e.g., see Reese et al. 2012). The second parameter, is 
the frequency at which detected modes of a star show the strongest 
amplitude. Brown et al. (1991) speculated that this characteristic 
frequency might scale with the atmospheric cut-off frequency for 
acoustic waves, Vac. Adopting an isothermal approximation to the 
full equation describing the cut-off frequency, one derives the scal¬ 
ing Vac Vmax (Brown et al. 1991, Kjeldsen & Bedding 

1995). Numerous studies have shown empirically that v^ax appears 
to follow this relation to reasonable approximation (e.g., Chaplin 
& Miglio 2013 and references therein). While some headway has 
been made on understanding the exact form of the scaling (e.g., 
Belkacem et al. 2011), a full theoretical justification remains elu¬ 
sive. Additions to the scaling, for example a dependence on the 
changing Mach number of the near-surface turbulent flow (Belka¬ 
cem et al. 2013), have also been proposed. 

Information encoded in Av and v^ax is commonly exploited by 
using the above scaling dependencies normalized to observed solar 
parameters or properties, i.e.. 
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Typical solar values for the seismic parameters (e.g., see Chaplin 
et al. 2014) are Avq = 135.1 yuHz and v^axo = 3090yuHz. These 
solar-calihrated scaling relations represent two equations in two 
unknowns when an estimate of is also available, allowing us 
to solve directly for M and R to give so-called “direct” estimates 
of the stellar properties. Moreover, each relation can be used inde¬ 
pendently to provide direct estimates of the mean stellar density p 
(Equation [TJ or the surface gravity g (Equation O. Alternatively, 
the relations may be utilized (together or separately) as part of a 
grid-based search code. Here, one searches amongst a grid of stel¬ 
lar evolutionary models to find those models whose predicted aster- 
oseismic or atmospheric parameters match the actual observed pa¬ 
rameters (at the level of the observational uncertainties). The solar- 
calibrated relations provide the means to translate model properties 
(R, M, Teff) to expected values for Av and v^.^, thereby allowing a 
comparison to be made with the observations. 

These global asteroseismic parameters and their associated 
scaling relations are now being employed in analyses of large sam¬ 
ples of solar-like oscillators to, for example, generate catalogues of 
asteroseismic stellar properties (e.g., see Huber et al. 2013, Chap¬ 
lin et al. 2014, Casagrande et al. 2014, Pinsonneault et al. 2014). 
Tests in the literature of the scaling relations for solar-type stars 
have, by and large, returned encouraging results. Studies have most 
commonly looked at data on estimated stellar radii, and include 
comparisons with very accurate properties from binaries, parallaxes 
and long baseline interferometry (e.g., Bruntt et al. 2010, Bedding 
2011, Miglio 2012, Huber et al. 2012, Silva Aguirre et al. 2012, 
White et al. 2013). Results have tested the combination of the two 
scaling relations to levels of around 4 per cent in inferred radii, and 
lOper cent in inferred masses (Chaplin et al. 2014). 

Results on red giants are more complicated. Eor example, He- 
core burning and H-shell burning giants with the same mass and 
radius can have a different Av, due to differences in the sound- 
speed profile in the outer layers (Miglio et al. 2012), which implies 
a different absolute scaling for Equation[T] Meanwhile, other stud¬ 
ies have looked at open clusters in the Kepler field, comparing re¬ 
sults on red giants inferred from asteroseismology and from turnoff 
eclipsing binaries (Brogaard et al. 2012, Sandquist et al. 2013). 

Our goal in this paper it to test, empirically, the accuracy of 
the classic scaling relation for oscillations seen in solar-type 
stars, i.e., cool main-sequence and sub-giant stars (we leave a study 
of giants to future work). There is actually very little in the liter¬ 
ature on the Vniax scaling alone, which partly reflects the difficulty 
of obtaining the data needed to test the one relation in isolation. 
A recent example made use of interferometric data on a few very 
bright Kepler targets; White et al. (2013) concluded that results on 
the E-type star Q Cyg may point to problems for the v^ax scaling in 
the hottest solar-type stars. 

Our basic approach is as follows. We use data on a sample 
of around 500 stars observed by the NASA Kepler Mission (the 
same sample as in Chaplin et al. 2014). Each star in our sample has 
a measured Vmax> which comes from analysis of the Kepler data. 
Throughout the rest of the paper we shall refer to these actual, mea¬ 
sured values as Vn,ax(data). Now, the classic scaling relation (Equa- 
tion|2j gives v^ax in terms of g and or, to be more specific, the 
combination g We can therefore test the scaling if we have 
independent measures of 
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to which the observed Vmax(data) may be compared. The Vmaxfgrid) 
are so-named because we adopt a grid-based search technique to 
estimate them, using as inputs the asteroseismic average large sep¬ 
arations, Av, of the stars along with photometric temperatures, T^g, 
derived using the Infrared Elux Method (IREM) and, where avail¬ 
able, metallicities [Fe/H] from spectroscopy. We search grids of 
stellar evolutionary models to find those models whose predicted 
|Av, Teif, [Fe/H]} match the actual observed inputs. Each model in 
the grid also has a computed v,„xx(grid). which comes from its M, 
R and T^g. The best-matching models will have the most likely val¬ 
ues of Vmax(grid). A suitable, likelihood-weighted average therefore 
provides an estimate of Vn,ax(grid) for every star (see Section |T2] |. 
which may then be compared directly to the observed Vn,ax(data). 
Any departures from a one-to-one correspondance of the values 
would point to problems with the classic scaling relation, and al¬ 
low us to quantify departures from the scaling. In sum, we leverage 
the potential of using asteroseismic results on a large number of 
solar-type stars to follow a statistical (ensemble) approach to the 
analysis and to thereby beat-down the errors. The approach is not 
dissimilar to that adopted by Morel et al. (2014) for part of their 
analysis of red giants observed by the CNES/ESA CoRoT Mission, 
which compared values of the surface gravity estimated using Av 
on the one hand and v^ax on the other. 

The layout of the rest of the paper is as follows. We begin 
in Section 12 with a discussion of the basic methodology adopted 
to test the Vn,ax scaling relation. Section [2 then introduces the real 
Kepler data, the artificial data we made to check our methodology 
(Section[3TT]l, and the grid-based search pipelines used to estimate 
the Vn,ax(grid) values (Section lT^ . Our results are presented in Sec- 
tion|4] beginning with a validation of the methodology using artifi¬ 
cial data (Section l4.lt . followed by results from the actual Kepler 
sample. We finish in Section[2with concluding remarks. 


2 METHOD 

We may in principle use some suitable grid-based results to test 
the Vn,ax scaling relation, albeit with caveats that we will discuss 
and address below. Let us suppose for the moment that grid-based 
searches using the set of inputs 

{Av, Teff, [Fe/H]j 

provide robust, unbiased estimates of the combination gT^g^, as 
calibrated to give Vn,ax(grid) defined by Equation [2 above. Assum¬ 
ing the temperatures and metallicities to be unbiased, at least to a 
level that will not influence significantly estimation of the combi¬ 
nation gT^g^, the fractional differences 

[Vmax(data)/Vmax(grid)] - 1, 

will provide a direct estimate of the bias in the Vmax scaling, i.e., 
the fractional amount by which the Vn,ax(data) values are over or 
underestimated relative to (g/go) {T^glTtSoY^^'^ v^axG- 

Crucial to the approach is the accuracy of Vn,ax(grid). First, we 
know to expect a small bias if the grid-based search pipelines em¬ 
ploy the Av scaling relation, as we now go on to explain. When the 
scaling relation is used in the grid-based searches, the fundamental 
properties of models in the grid are used as inputs to Equation [T] 
to yield model estimates of Av for comparison with the observed 
separations. Alternatively, one may circumvent use of the scaling 
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by computing for each model a set of theoretical oscillation fre¬ 
quencies (e.g., radial-mode frequencies spanning the same orders 
as those observed in the real data), from which one may then esti¬ 
mate the required Ay from a suitable fit to those frequencies. 

It is now well known that predictions made by the calibrated 
scaling-relation (Equation [TJ have small, systematic differences 
with respect to predictions from model-computed frequencies (e.g., 
see Ulrich 1986, White et al. 2011; Mosser et al. 2013). For solar- 
type stars, these differences can be up to - 2 percent in size, and 
become more pronounced at effective temperatures progressively 
further away from ^ 5700 K (e.g., see Figures 5 and 6 of White 
et al. 2011; and figures in Chaplin et al. 2014). Here, we have em¬ 
ployed grid-based search pipelines that can run with or without the 
Ay scaling relation. We use one pipeline that may he run either 
using individual model-calculated frequencies or the Ay scaling re¬ 
lation; and, for comparison, two other pipelines that used the Ay 
scaling only. 

Second, we must consider the impact of the poor modelling 
of the near-surface layers of stars. In the case of the Sun it is now 
well established that this gives rise to a frequency-dependent off¬ 
set between observed and model-calculated oscillation frequencies 
(e.g., see Chaplin & Miglio 2013, and references therein). This so- 
called “surface term” increases in magnitude with increasing over¬ 
tone number, n. The amount by which Ay is affected will depend 
on the variation of the surface term with n. Tests of the grid-based 
method (Basu et al. 2010) indicate that the impact of the solar sur¬ 
face term on Ay - which decreases the observed solar Av by just un¬ 
der 1 per cent compared to model predictions - leads only to very 
small errors in the inferred solar properties, certainly well within 
the observational uncertainties associated with the Kepler data used 
in this paper. Stellar surface terms would need to be substantially 
larger than the solar term to produce significant bias in our results. 
However, the nature of the term in other stars remains rather poorly 
understood, and so this caveat should be borne in mind. 

Our tests with artificial data do nevertheless provide some in¬ 
sights on the sensitivity of the results to such offsets. The artificial 
data come from models computed using a different stellar evolu¬ 
tionary code and different input physics to those of the grids to 
which the grid-based pipelines are coupled. This can give rise to 
“surface-term” like offsets between models in the artihcial sample 
in the grids that share the same fundamental properties (in partic¬ 
ular from differences in boundary conditions, matching to model 
atmospheres etc.). 

3 DATA AND GRID PIPELINES 
3.1 Real and artificial data 

The observational data for our study come from Chaplin et al. 
(2014). This study produced an asteroseismic catalogue from an 
extensive grid-based analysis of more than 500 solar-like oscilla¬ 
tors, which were observed by Kepler as part of an asteroseismic 
survey that was conducted over the first 10 months of science oper¬ 
ations. Stellar properties were estimated from a grid-based analysis 
using the global asteroseismic parameters Av and v^ax together with 
complementary photometric and spectroscopic data as the inputs. 
Homogeneous sets of effective temperatures were available for 
the full sample of stars, courtesy of complementary ground-based 
photometry. A homogeneous set of spectroscopic parameters (Tefr 
and [Fe/H]) was also available from Bruntt et al. (2012), but only 
for a subset of 87 stars in the sample 


Here, we make use of the global asteroseismic parameters 
and the complementary data to perform the new grid-based anal¬ 
ysis needed to test the y^ax scaling relation. The Chaplin et al. 
(2014) sample is dominated by cool main-sequence and sub-giant 
stars but does contain a small fraction of stars at the base of the 
red-giant branch that were serendipitously observed as part of the 
short-cadence asteroseismic survey. We have removed the more 
evolved stars from the sample (which will be the subject of a sepa¬ 
rate study). The selected sample contains 426 solar-type stars. 

As noted above, complementary photometry was available on 
the entire sample. This allowed us to perform a new, homogenous 
grid-based analysis on all the selected solar-type stars, but at the 
cost of not having robust, well-constrained estimates of [Fe/H] for 
each star since the complementary photometry available to us in 
the Kepler Input Catalogue (KIC; see Brown et al. 2011) does 
not provide strong constraints on metallicity. Just like Chaplin et 
al. (2014), we therefore adopted an average [Fe/H] value as input 
for every star when we analysed the 426-star sample. We actually 
tried two different values; One set of results came from using the 
average of the 87 metallicities measured by Bruntt et al. (2012), 
i.e., [Fe/H] = -0.05 dex; while the other results came from us¬ 
ing [Fe/H] = -0.20 dex for all stars (e.g., see Silva Aguirre et al. 
2011), the value adopted in Chaplin et al. (2014). In both cases we 
adopted large input uncertainties of ±0.3 dex. 

In spite of the weak constraints on [Fe/H], we still obtained 
more precise results from the larger 426-star sample having com¬ 
plementary photometric data than we did from the smaller sample 
with complementary spectroscopic data because the larger sample 
size compensated for the inferior precision in [Fe/H]. Results ob¬ 
tained were similar, and hence in what follows we present detailed 
results from the photometric sample. This sample also provided 
much better coverage in the domain where Tf.g > 6000 K. 

Finally with regards to the input data, we note that the pho¬ 
tometric temperatures were the Infrared Flux Method (IRFM) es¬ 
timates from Chaplin et al. (2014), which were calculated us¬ 
ing multi-band photometry in the Two Micron All Sky Survey 
(2MASS; see Skrutskie et al. 2006) JHK bands, and in the Sloan 
Digital Sky Survey (SDSS) griz bands (both available in the KIC). 

Panel (a) of Fig. [T] plots the locations of the sample of se¬ 
lected Kepler solar-type stars on a Hertzsprung-Russell diagram. 
The lines mark evolutionary tracks computed by the Modules for 
Experiments in Stellar Astrophysics (MESA) code (Paxton et al. 
2011; see below), for models having solar composition and masses 
ranging from 0.8 to 1.5 Mq (in steps of 0.1 Mq). Panel (b) of the 
same figure shows the measured Vn,ax(data) of the Kepler sample, 
from analysis of the oscillation spectra of the stars, indicating the 
parameter range tested by the analysis. 

We also used artificial data to test and validate our method¬ 
ology. The artificial sample of stars was comprised of models 
drawn from an evolutionary grid, computed using MESA. The grid 
spanned the range 0.8 < M/Mq <1.5 and -0.8 < [Fe/H] < 0.8 
(both in steps of 0.1), with tracks computed from the pre-main se¬ 
quence to the base of the red-giant branch. 

The Grevesse & Noels (1993) value of Zq/Xq = 0.0245 was 
used to translate between model values of [Fe/H] and Z/X. The 
OPAL equation of state (Rogers & Nayfonov 2002) and OPAL 
opacities (Iglesias & Rogers 1996) were used, augmented by low- 
temperature opacities from Eerguson et al. (2005). Nuclear reac¬ 
tion rates were from NACRE (Angulo et al. 1999), with updates 
for the ''*N(p, y)'^0 (Imbriani et al. 2004, 2005) reaction. Con¬ 
vection was treated according to mixing-length theory, using the 
solar-calibrated mixing length parameter. Diffusion and effects of 
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Figure 1. Panel (a): Hertzsprung-Russell diagram of the asteroseismic sample of Kepler solar-type stars. Panel (b): measured values, Vmax(data), from analysis 
of the oscillation spectra of the real stars, also as a function of effective temperature. Panel (c): Hertzsprung-Russell diagram of the sample of artificial targets, 
computed from MESA models. The lines in the panels show evolutionary tracks computed by the MESA code (see text), for models having solar composition 
and masses ranging from 0.8 to 1.5 Mq (in steps of 0.1 Mq). 


rotational mixing were not included. The primordial Helium abun¬ 
dance was fixed to Tp = 0.2484 (Cyburt et al. 2003), and the helium 
enrichment set to AT/AZ = 2 (e.g., Chiosi & Matteucci 1982, Ca- 
rigi & Peimbert 2008). We refer the reader to Paxton et al. (2011) 
for further details. 

Artificial stars were drawn from the grid by seeking a “best 
matching” model for each of the solar-type stars in the Kepler sam¬ 
ple (via diX^ minimization). We found that a selection based solely 
on a comparison of the observed {Kepler) and model (grid) values 
of Tj-ff and Av (the latter using the calibrated scaling relation) was 
sufficient to produce an artificial sample that had distributions in 
each of the fundamental properties that were a reasonable match to 
those for the real sample. Panel (c) of Fig.[T]marks the locations of 
the selected models in our artificial sample of stars (using the pris¬ 
tine model parameters). Owing to the limited resolution of the grid, 
we found that some real stars had the same best-matching artificial 
model. This meant that our final sample of selected artificial stars 
was comprised of 306 unique models. 

Next, we computed adiabatic oscillation frequencies for each 
of the 306 selected models, using the GYRE (Townsend & Teitler 


2013) stellar oscillations code. The Av of each artificial star was 
then given by the best-fitting gradient of a linear fit to the radial 
order, n, of the five / = 0 frequencies centred on the estimated 
(see below) of the model. This approach gives values that are rep¬ 
resentative of the average values extracted from the real data. We 
also tested the impact of calculating Ay using a different number of 
orders (to reflect the varying data quality and S/N levels in the real 
sample of stars) and from adopting a weighted fit of the frequencies 
(to reflect the impact of the changing S/N in any given observed 
spectrum). Neither change had a significant impact on our results. 

We computed three different sets of y,„ax for the artificial sam¬ 
ple. One set was computed assuming perfect adherence to the solar- 
calibrated scaling relation, i.e., by using Equation|2] We made two 
other sets by applying a temperature-dependent fractional offset to 
the solar-calibrated scaling relation, i.e., a computation that took 
the form 

v.„ = r(reff)(— j|-J y_e, (4) 

with 'F(7’e(f) being the fractional temperature-dependent offset. We 
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applied a linear and quadratic offset, respectively, both of which 
are shown in panel (a) of Fig. |4] (in the figure that also shows re¬ 
sults from the data; see below for further discussion). Use of the 
first set of “perfect” values allowed us to test the impact of biases 
not associated with the v^ax scaling relation; while the other sets al¬ 
lowed us to test whether we could recover information on a known 
bias in the v^ax scaling. 

Finally, with each artificial star then having calculated values 
of Av and v^ax, and a model it remained to add noise and to as¬ 
sign uncertainties to those parameters for input to the grid pipelines. 
This meant that in our analysis we would treat the artificial data in 
exactly the same way as the real data. Flere, we simply used the 
relevant parameter uncertainties of the real Kepler star to which 
each artificial star was associated. To make a given realization of 
the artificial datasets, we added Gaussian noise to each pristine in¬ 
put parameter, multiplied by the relevant parameter uncertainty, to 
give Vn,ax(data), Av(data), and TeffCdata) for each artificial star. 

3.2 Grid pipelines 

We utilized three different grid-based pipeline codes to return esti¬ 
mates of the parameter Vn,ax(grid) for all stars in the real and artifi¬ 
cial samples: 

- the Bellaterra Stellar Properties Pipeline (BeSPP) (Serenelli et 
al. 2013, extended for asteroseismic modelling); 

- the Yale-Birmingham (YB) (Basu et al. 2010, 2012, Gai et al. 
2011); and 

- PARAM (da Silva et al. 2006; Miglio et al. 2013; Rodrigues 
et al. 2014); 

The BeSPP pipeline was run with a grid comprised of models con¬ 
structed with the GARSTEC code (Weiss & Schlattl 2008). The 
parameters of the grid are described in Silva Aguirre et al. (2012). 
BeSPP was run in two different modes of operation: one where 
grid-model estimates of Av were computed using adiabatic oscil¬ 
lation frequencies (frequency mode); and one where the estimates 
were instead computed using the solar-calibrated scaling relation 
(Equation[TJ scaling-relation mode). The other two pipelines were 
run only in the latter, scaling-relation mode. PARAM was run us¬ 
ing a grid comprising models made by the Padova group (Marigo 
et al. 2008). Eurther details may be found in Miglio et al. (2013). 
The YB pipeline returned results using five different sets of stel¬ 
lar models: grids computed by the Dartmouth (Dotter et al. 2008) 
and Padova groups (Bressan et al. 2012); the set of YY isochrones 
(Demarque et al. 2004); a grid comprised of the BASTI models of 
Pietrinfemi et al. (2004), computed for use in asteroseismic studies 
(see Silva Aguirre et al 2013); and, finally, a grid constructed us¬ 
ing the YREC Code (Demarque et al. 2008), which is described by 
Basu et al. (2012). This grid has been used in other papers, and we 
retain the YREC2 name here. 

We also report results from YB which are labelled ALL. This 
set of results was generated by combining the YB analysis over all 
five grids to compute what are in essence averages from a compos¬ 
ite “super distribution”. This is a new addition to the YB code that 
has not previously been documented in the literature, and so we 
provide further details here. 

The YB pipeline determines the properties of a star using the 
given observational input (central) parameter set. A key step in the 
method is to generate 10,000 input parameter sets by adding dif¬ 
ferent random realisations of Gaussian noise (commensurate with 
the input uncertainties) to the actual (central) observational input 
parameter set. For each realisation, we find all models in a grid 


within 3cr of the input uncertainties, and use these models to define 
a likelihood function (e.g., see Gai et al. 2011, Basu et al. 2012 
for details). The estimated property, e.g., Vn,ax(grid), is then the 
likelihood-weighted average of the property of the selected models. 
The 10001 values of any given property estimated from the central 
value and the 10000 realisations, form the probability distribution 
function for that parameter. In the YB pipeline we adopt the median 
of the distribution as the estimated value of the property, and we use 
Icr limits from the median as a measure of the uncertainties. The 
ALL results were obtained by constructing a consolidated probabil¬ 
ity distribution for a given star by adding together the five different 
distribution functions obtained using the five different grids. Then 
we determined the median of this consolidated or “super” distribu¬ 
tion function, to determine the average Vmaxfgrid). 

All three pipelines were employed in the grid-based analysis 
described in Chaplin et al. (2014), where summary details of the 
physics employed in the grids may also be found. 

4 RESULTS 

4.1 Results from artificial data 

We begin with results from analysing the first set of artificial data. 
Recall that in this case v^axfdata) for each artificial star - essentially 
the proxy for what would be the observed v^ax of each real star 
- was computed assuming strict adherence to the solar-calibrated 
scaling relation. The top two panels of Fig. |2] show results from 
the BeSPP pipeline, run in both frequency mode [panel (a)] and 
scaling-relation mode [panel (b)], on a representative noise re¬ 
alization of the artificial set. Both panels show fractional differ¬ 
ences Vmax(data)/vmax(grid) - 1, i.e., fractional differences between 
the simulated measurements, Vn,ax(data), and the grid-based esti¬ 
mates, Vn,ax(grid). Here, grid-based estimates were calculated as¬ 
suming an input [Fe/H] pegged to the Bruntt et al. average, i.e., 
[Fe/H] = -0.05 dex. Results are shown both for the individual arti¬ 
ficial stars (symbols) and for averages computed over 130-K bins in 
Teif (lines). The error bars mark uncertainties on each average. The 
scatter in the fractional differences is entirely consistent with the 
formal uncertainties on the differences (which are propagated from 
the individual formal uncertainties on Vn,ax(data) and Vn,ax(grid)). 

The results on the whole show the expected trends. The 
scaling-relation-mode results in panel (b) are not flat. The differ¬ 
ences show an upward trend with increasing T eff, which is due to 
the known offsets in the Av scaling. This tells us that if there were 
no Vn,ax bias in the real data we should expect to see trends like 
these when using a scaling-relation-based grid pipeline. 

With regards to the frequency-mode results in the panel (a) of 
Fig.H as expected the upward trend from the Av scaling is absent, 
and we see a flat trend in the comparison at the level of precision 
of the data. The frequency-mode results therefore allow us to con¬ 
clude, correctly, that the artificial data follow a like scaling. 

That said, the differences do show a small absolute offset, albeit 
one that does not change significantly with TcB- The absolute offset 
has more than one contribution. 

First, there is a contribution due to the uncertainty in [Fe/H] 
we adopted for the input data (to mimic that for the real 426-star 
sample). Panel (c) of Fig. shows results from BeSPP run in fre¬ 
quency mode, but with [Fe/H] = -0.20 dex now used as input for 
all stars. This change to the input [Fe/H] produces a shift in the 
absolute offset of just under 1 per cent. Second, there will be a con¬ 
tribution to the absolute offset arising from differences between the 


6 Coelho et al. 


(a) (b) 




(c) 



Figure 2. Results from BeSPP grid pipeline, for artificial data following perfect adherence to the Vmax scaling relation. Panel (a) shows results from when 
BeSPP is run in frequency mode; while panel (h) shows results when the pipeline is run in scaling-relation mode (see text). Plotted are fractional differences 
rmax(data)/Vn,ax(grid) — 1, with grid-based estimates calculated assuming an input [Fe/H] of—0.05 dex for every artificial star. Results are plotted for individual 
artificial stars (symbols) and for averages computed over 130-K bins in T^g (lines). The eiTor bars mark uncertainties on each average. Panel (c): Frequency¬ 
mode results from analysing the same artificial data, but now with an input [Fe/H] of —0.20 dex. 


models on which the artificial stars are based and those used in the 
BeSPP grid. These differences can give rise to what looks like a 
“surface term” offset (see discussion in Sectionabove). We con¬ 
clude that we should expect there to be some uncertainty over the 
absolute calibration of the relation for results on the real Kepler 
data. We shall come back in Section |4^ to attempt an estimate of 
the relative contributions of the above effects to the uncertainty in 
the absolute calibration. 

Fig. H also includes results from the YB and PARAM 
pipelines, which were coupled to a variety of grids with [Fe/H] = 
-0.05 dex used as input. The YB ALL results are plotted in black, 
the YB results from individual grids in grey, and the PARAM re¬ 
sults in blue. Recall that these pipelines were run only in scaling- 
relation mode. To aid the clarity of the plots we present just the 
Feff -binned averages. Results from the other pipelines follow the 
BeSPP scaling-relation results from the top panels of Fig.[2 which 
we plot again here (in red) for direct comparison. The shapes of 
the scaling-relation trends for each grid-pipeline combination are 
similar - showing an upward trend with increasing T^g, again due 


to the known offsets in the Av scaling - with an extreme spread of 
approximately 3.5 per cent in the fractional differences. 

Fig-Elpresents results for the artificial sets which have known 
biases in Panel (a) shows the linear and quadratic biases, 
'F^(Fcff) (Equation |4l(, which we imposed on the artificial data. The 
other panels show results of testing the biased artificial data, with 
fractional differences plotted for each biased set (see figure leg¬ 
ends). Also plotted for reference (black lines) are the recovered 
trends from the unbiased data (see Figs.j^andO. The panels show 
returned fractional differences from: BeSPP in frequency mode (b); 
BeSPP run in scaling mode (c); YB ALL (d); and PARAM (e). 

The results are again encouraging. Evidence of the injected 
bias is clearly present in the results. We also see the offset between 
the two bias trends, although we would not be able to tell the differ¬ 
ence between the shapes of the trends. Nevertheless, we are able to 
conclude that for the level of bias tested here, it would be possible 
to discriminate between the no-bias and bias cases. 
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Figure 3. Fractional differences Vmax(data)/vmax (grid) — 1 returned by all three grid pipelines, for aitiffcial data following perfect adherence to the solai*- 
calibrated v^ax scaling relation (Equation^. Different lines show 130-K averages in Teff for different pipeline-grid combinations (see annotation), all of which 
used an input [Fe/H] of —0.05 dex. 


4.2 Results from Kepler data 

Fig. a shows results from analysing the selected sample of Kepler 
solar-type stars, which all have complementary photometric data. 
The fractional differences Vn,ax(data)/Vn,ax(grid) - 1 are plotted for 
the different pipelines, including results obtained using BeSPP in 
both modes of operation. The top four panels show results with 
[Fe/H] = -0.05 dex used as input, as returned by: BeSPP in fre¬ 
quency mode (a); BeSPP run in scaling mode (b); YB ALL (c); 
and PARAM (d). 

The trends observed here bear a striking resemblance to those 
given by the results on the artificial data that follow the classic v^ax 
scaling relation. The trend in the BeSPP frequency-mode Kepler 
results on the real data is flat in Feffi and consistent with a 
like scaling at the level of precision of the data. 

The scatter in the fractional differences of the real data in 
Fig.S is entirely consistent with statistical scatter, given the for¬ 
mal uncertainties. The histograms in Fig. show the normalized 
distributions of fractional differences for fhe real (red) and arti¬ 
ficial no-bias (blue) BeSPP frequency-mode results. Here, each 
t'max(data)/Vn,ax(grid) - 1 was subtracted from the respective mean 
trend-line of its sample, and then normalized by its formal uncer¬ 
tainty, that uncertainty having been propagated from the individual 
formal uncertainties on Vm^xCdata) and Vn,ax(grid). The lines show 
the cumulative distributions of the histogram data. The real and ar¬ 
tificial data histograms follow one another very closely. 

We may quantify further the adherence of the BeSPP 
frequency-mode results to a like scaling, i.e., a flat trend. 

If we fit a simple linear model to all the data, the best-fitting gradi¬ 
ent implies a linear bias of 0.16 ± 0.10 per cent per 100 K. We may 
therefore rule out departures from a gT^f^ scaling at the level of 
=^1.5 percent over the full AT =; 1560 K range tested in Feff. 

As noted in Section[T] the scaling relations (Equations[T]and[2j 
may be manipulated to give expressions for the stellar radius, R, and 
mass, M. The dependencies of the resulting expressions imply that 
any bias in Vn,ax propagates to give bias in the inferred properties of 



and 

lSM\ ^ / dVmax \ 

' M / i V„ax j ■ 

Departures from the classic scaling of =: 1.5 percent across the full 
range (see above) therefore fix 1-cr upper limits on changes in any 
bias in inferred properties with changing Teff of SR/R <1.5 percent 
and SMjM < 4.5 percent. Since the v^ax scaling relation implies 

/^\ ^ / dVmax \ 

\ g / \ Vmax / ’ 

bias in surface gravities inferred from use of the relation will be 
similar to that for the radii above (which corresponds to a bias in 
log g of =: 0.006 dex). 

We also fitted data in overlapping ATefr = 500 ^-wide bins, 
moving systematically through the full range to test the temperature 
dependence of any departures from the scaling. Throughout most of 
the range we obtain 1-cr uncertainties on the best-fitting gradients 
of - 0.4 per cent per 100 K. The limits increase to 0.5 per cent per 
100 K at temperatures above about 6200 K, rising to ^ 1 per cent 
per 100 K in the highest range we fitted, which started at 6450 K. 

Our results neverthless imply some uncertainty over the ab¬ 
solute calibration of the v^ax scaling. Panel (e) in Fig. |5] shows the 
impact on the BeSPP frequency-mode results of using a different 
input [Fe/H] for each star, here -0.20 dex instead of the value of 
-0.05 dex used for the top left-hand panel, lust like the artificial 
data in Fig. [2] we see a systemetic shift in the absolute otfset, from 
just over 3 per cent for an input [Fe/H] of -0.05 dex to around 4 per 
cent for an input [Fe/H] of -0.20 dex. Assuming that the overall 
shift results from the quadratic addition of uncorrelated effects, this 
implies a systematic shift due to the uncertainty in [Fe/H] of around 
2.5 per cent. It is worth noting that if we analyse the smaller sample 
of real stars having spectroscopic data from Bruntt et al. (2012) - 
which have much more tightly constrained input [Fe/H] - we get 
the same 3 per cent offset as the -0.05 dex full-sample case. 

Accounting for the above still leaves us with an uncertainty in 
the absolute calibration of around 3 per cent (again, assuming that 
the contributions are uncorrelated). We recall from discussions in 
previous sections that we would expect this to reflect differences 




































(data)/y„„(grid) - 1 y„^(data)/y„^(grid) 
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(a) 



(b) 


(c) 




(d) 




Figure 4. Results from all grid pipelines, using an input [Fe/H] of —0.05 dex, for artificial datasets having a Te^-dependent bias imposed on the solar-calibrated 
Vmax scaling relation. Panel (a): Linear and quadratic biases, (Equation |4), imposed on the artificial data. Other panels: results of testing the biased 

artificial data. Plotted are the fractional dilferences Vmax(data)/Vmax(grid) — 1 for each biased set, returned by: BeSPP in frequency mode (b); BeSPP nm in 
scaling mode (c); YB ALL (d); and PARAM (e). Each panel also shows for comparison the result from the unbiased data (see Figs.[2]an dE). 
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Figure 5. Results from analysis of the real Kepler sample of solar-type stars, showing the fractional differences Vmax (data)/v^ax(grid) - 1 returned by the 
different pipelines. The top four panels show results with [Fe/H] = —0.05 dex used as input, as returned by: BeSPP in frequency mode (a); BeSPP run in 
scaling mode (b); YB ALL (c); and PARAM (d). Panel (e): BeSPP frequency-mode results form using an input [Fe/H] of —0.20 dex for every star. 
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Figure 6. Bar plots (associated to left-hand ordinate scale): histograms of the normalized distribution of fractional differences Vmax(data)/Vmax(grid) “ 1 for 
the artificial no-bias (blue) and real (red) BeSPP frequency-mode results, which were plotted in Figs.|2]and|^respectively. Each fractional difference was 
subtracted from the respective mean trend-line of its sample, and then normalized by its formal uncertainty, that uncertainty having been propagated from the 
individual formal uncertainties on Vn,ax(data) and Vn,ax(grid)- Lines (associated to right-hand ordinate scale): the cumulative distributions of the histogram data. 


between the physics and structures of the real stars (including near¬ 
surface effects) and those of the grid models. 

An analysis performed on ATeff = 500 A'-wide bins finds no 
evidence for any significant deviation of the offset with The 
results imply that if there are any absolute errors in the calibration, 
all inferred properties will be biased by the same fractional amount. 
For the solar-calibrated Equation the calibration is provided by 
the value VmaxG = 3090yuHz. Guided by the results above, a 3 per 
cent uncertainty in the overall calibration translates to an error in 
the the calibrating frequency of around lOOyuHz. 


5 CONCLUSION 

Asteroseismic scaling relations for solar-like oscillators are becom¬ 
ing increasingly important as they receive more widespread use in 
building large catalogues of stellar properties. 

We have tested the scaling relation for the global asteroseis¬ 
mic parameter v^ax, the frequency at which a solar-like oscilla¬ 
tor presents its strongest observed pulsation amplitude. The classic 
scaling relation assumes that Vmax scales with surface gravity and 
effective temperature according to We have tested how well 

the detected oscillations in a large sample of solar-type stars ob¬ 
served by the NASA Kepler Mission adhere to this relation by com¬ 
paring the observed Vmax of the stars with independent estimates of 
the combination g 

Our results rule out departures from the classic v^ax scaling at 
the level of 1.5 percent over the full =; 1560 K range in that 
we tested. There is some uncertainty over the absolute calibration 
of the scaling, but any variation with is evidently small, with 
limits similar to those above. 
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